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Abstract 

We report a measurement of the inclusive electron energy spectrum for charmed semileptonic 
decays of B mesons in a 140 fb -1 data sample collected on the T(4S) resonance, with the Belle 
detector at the KEKB asymmetric energy e + e~ collider. We determine the first, second and third 
moments of the electron energy spectrum for threshold values of the electron energy between 0.4 
and 1.5 GeV. 

PACS numbers: 12.15.Hh, 11.30.cr, 13.25.Hw 
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INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa matrix element \V c b\ can be extracted from the in- 
clusive branching fraction for charmed semileptonic I?-meson decays B(B -> X c £i/)0, Q. 
Several studies have shown that the spectator model decay rate is the leading term in a 
well-defined expansion controlled by the parameter Aqcd/"V Non-perturbative corrections 
to this leading approximation arise only at order 1/m^ and above. The key issue in this 
approach is the ability to separate non-perturbative corrections, that can be expressed as a 
series in powers of l/m&, and perturbative corrections, expressed in powers of a s . 

The coefficients of the l/m& power terms are expectation values of operators that in- 
clude non-perturbative physics. Different expansions exist, reflecting a difference in the ap- 
proach used to handle the energy scale /x which separates long-distance from short- distance 
physics jH- 

The shape of the lepton spectrum provides constraints on the heavy quark expansion 
based on local Operator Product Expansion (OPE) jiJ], which calculates properties of the 
B — > X c £u transitions. So far, measurements of the electron energy distribution have been 
made by DELPHI, CLEO, BABAR and Belle collaborations 0, IS In tnis note we 
report a measurement of the first, second and third moment of the electron energy spectrum 
with a minimum electron momentum cut ranging between 0.4 and 1.5 GeV/c in the B meson 
rest frame. 

The data used in this analysis were collected with the Belle detector at the KEKB jj| 
asymmetric energy e + e~ collider. The Belle [HI detector is a large-solid-angle magnetic 
spectrometer that consists of a three- layer silicon vertex detector (SVD), a 50-layer central 
drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel-like 
arrangement of time-of-flight scintillation counters (TOF), and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals (ECL) located inside a super-conducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return located outside of the coil is instru- 
mented to detect K\ mesons and to identify muons (KLM). 

The present results are based on a 140fb _1 data sample collected at the T(45") resonance 
(on- resonance) , which contains 1.5 x 10 8 BB pairs. An additional 15 fb — 1 data sample 
taken at 60MeV below the T(45") resonance (off-resonance) is used to perform background 
subtraction from the e + e~ — > qq process. Events are selected by fully reconstructing one of 
the B mesons, produced in pairs from T(4S) decays. 

We used a fully simulated generic Monte Carlo sample equivalent to 2.4 times the on- 
resonance integrated luminosity. Simulated Monte Carlo events are generated with the 
EVTGEN event generator, and full detector simulation based on GEANT is applied 
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EVENT SELECTION 

After selecting hadronic events 0], we fully reconstruct the decay of a B meson on one 
side (tag-side), in the decay modes B — > D^*'n + , D^*> p + , D^af , yielding a high purity B 
meson sample. The following sub-decay modes are reconstructed: 

. D*° — > D°7T°, Z)°7, 

• D*~ -> D°n + ,D^n°, 

• D° — > Kti, Ktiti , Ktititi, Ksim, Ksii and 



1 



• D — > Ktxtx, Kstt. 

For each selected event, we calculate the beam-constrained mass Mb c and energy differ- 
ence AE: 

M hc = ^/(^eanJ 2 " iP%) 2 , AE = E* B - E^ (1) 

where E£ eam , p* B and E* B are the beam energy, the B momentum and the B energy in the 
centre of mass frame, respectively. Events with M^ c > 5.27 GeV/c 2 and —0.06 GeV < AE < 
0.08 GeV are considered to be signal candidates. The number of events, after continuum 
subtraction, in the signal region are 63155 ± 931 and 40032 ± 475, B + and B° candidates, 
respectively [l3j |. 

ELECTRON SELECTION 

We search for electrons produced in semileptonic B decays on the "non-tag" side. Particle 
tracks are selected if they originate from near the interaction vertex. In addition, we measure 
tracks which pass through the barrel region of the detector, corresponding to an angular 
acceptance of 35° < 6*i a b < 125°, where 9^ is the polar angle of the track relative to the z 
axis (opposite to the positron beam line). 

Tracks that pass the above selection criteria and are not used in the reconstruction 
of the tag-side B meson are considered as electron candidates. Electron identification is 
based on a combination of ionisation dE/dx measurements in the CDC, the response of 
the ACC, the shower shape in the ECL, and the ratio of energy deposited in the ECL 
to the momentum measured by the tracking system (E/p) |3~ol | - The electron momentum 
spectrum is corrected by a momentum dependent electron detection efficiency, which includes 
the detector acceptance as well as the selection efficiency The momentum of the selected 
electrons is calculated in the B meson rest frame {pl B ), exploiting the knowledge of the 
momentum of the fully reconstructed B. We require p* e B > 0.4 GeV/c. The electron yield 
after these cuts is given in Table |H with a signal purity of 69.6% (50.0%) for B + (B°) tags. 

Electrons suffer a loss of energy due to bremsstrahlung radiation in the detector material 
in front of the calorimeter. This biases the electron energy distribution. Hence, we partially 
correct for this energy loss by recovering some of the emitted photons to restore the original 
electron energy. The photon candidate is combined with the electron if the photon energy 
is below 1 GeV and the angle between the photon and the electron is less than 0.05 radians. 

In B + decays, prompt semileptonic decays (b — > xlv) of the non-tag side B mesons are 
separated from cascade charm decays (b — > c — > ytv), based on the correlation between the 
flavor of the tagging B and the electron charge. In B° decays, part of the sample mixes, 
which flips the correlation to the B°. Thus in the B° sample we do not cut on the electron 
charge — B flavor correlation. 

BACKGROUND SUBTRACTION 

The reconstructed electron energy spectrum is contaminated by background processes, 
which should be evaluated and subtracted from the distribution before the extraction of the 
moments. The residual background is from: 

• continuum background, 
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• combinatorial background, 

• background from secondary decays, 

• J/x/j,ip(2S), Dalitz decays and photon conversions, 

• fake electrons. 

The shape of the continuum background is derived from off-resonance data, and is nor- 
malised using the off- to on-resonance luminosity ratio and cross section difference. To 
account for low statistics in the off-resonance data we fit an exponential to the distribution, 
before rescaling. The shape of the combinatorial background is derived from generic BB 
Monte Carlo events where either reconstruction or flavor assignment of the tagged B meson 
is not carried out correctly. The yield of this background is normalised to the on resonance 
data Mb c side band (Mb c < 5.25 GeV/c 2 ), after the continuum background subtraction. 
We also correct for cases where the fully reconstructed B is correctly tagged, but the elec- 
tron candidate either is not from a B decay (secondary) or is a mis-identified hadron. These 
background sources are irreducible; to estimate their magnitude we normalise the BB Monte 
Carlo using a fit to the electron momentum distribution after continuum and combinatorial 
background subtraction. 

The background distribution from B — > D^-*' — > e decays are scaled using the latest pub- 
lished branching fractions jT3j- Contributions from J/ip,ip(2S) decays, photon conversions, 
and Dalitz decays are small after the track and electron selection cuts. The surviving back- 
grounds are estimated using Monte Carlo simulation and subtracted along with the major 
secondary backgrounds. The normalisation for the Monte Carlo yield of secondary and fake 
leptons is obtained from data by fitting to the lepton momentum distribution in the range 
0.3 GeV/c < pf < 2.4 GeV/c. Figure □ shows the raw electron momentum spectrum with 
all background contributions overlaid. Table |1] summarises the number of detected electrons 
and the contributions from these backgrounds. 

TABLE I: Electron yields for p* e B > 0.4 GeV/c. The errors are statistical only. 



B candidate 




B+ 




B° 


On Resonance Data 


6573 


± 


81 


5564 


± 75 


Scaled Off Resonance 


258 


± 


48 


218 


± 45 


Combinatorial Background 


1394 


± 


38 


765 


± 28 


Secondary 


680 


± 


26 


1915 


± 44 


Background Subtracted 


4241 


± 


105 


2666 


± 102 



THE ELECTRON ENERGY SPECTRUM AND THE MOMENTS 

The electron energy spectrum is generated via Monte Carlo simulation of B — > X c ev 
decays with the EVTGEN event generator 1 1 111 . The spectrum from B — > X c eu is modelled 
using four components: X c = D (ISGW2 D* (HQET ^H), higher resonance charm 

meson states D**(ISGW2) and non-resonant D^n (Goity and Roberts 0]). To account for 
the most recent theoretical and experimental results, we reweight the D and D* components 
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FIG. 1: Breakdown of the backgrounds in the electron momentum spectra for B + , and B° elec- 
trons. The Monte Carlo sample does not include b — > u transitions. The errors shown are statistical 
only. 



in pi to the spectra generated with current (world) average form factors [17J • In addition, 
the branching fractions for the D and D* components are corrected according to current 
(world) average values [l7j |. Electrons that come from the b —> u transition are subtracted 
from the unfolded electron energy spectrum, defined later. We model the electron energy 
spectrum from B — > X u lu transitions using the De Fazio and Neubert prescription ^8[. The 
6-quark motion parameters are derived in Ref. 0|. We scale according to the B — > X u lv 
branching fraction in Ref. (l7| . 

To measure the first, second and third electron moments we need to determine the true 
electron energy spectrum in the B meson rest frame, E* B . The background subtracted mo- 
mentum spectrum is distorted by various detector effects. The true electron energy spectrum 
is extracted by performing an unfolding procedure based on the Singular Value Decompo- 
sition (SVD) algorithm |20j]. The unfolded spectrum is corrected for QED radiative effects 



using the PHOTOS algorithm |21|, as the OPE does not have an 0(a) QED correction. 
The unfolded electron energy spectrum is shown in Figure 121 

We measure the central moments of the electron energy spectrum. The first moment is 
defined to be M[ = (Ej) and subsequent central moments are determined about the first 
moment, = ((Ei — (Ei)) n ), where I is the electron energy cutoff and n — 2,3. We 
measure the first three moments with six electron energy threshold cuts (/ = E cnt = 0.4, 
0.6, 0.8, 1.0, 1.2 and 1.5 GeV) combining the spectra from B + and B° semileptonic decays. 
Table |H] provides the final measurements of the moments. Figure H3 shows the moments of 
the B° and B + subsamples as a function of E cnt , as well as the B° and B + combined average. 
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E* B (GeV) 



FIG. 2: Unfolded electron energy distribution in the B meson rest frame, combining contributions 
from B° and B + decays, and corrected for QED radiative effects. The errors shown are statistical. 

SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties in the moments stem from event selection, electron identi- 
fication, background estimation and model dependence. 

We determine the electron identification uncertainty from varying the electron identifi- 
cation constraints. In addition, we calculate a systematic uncertainty associated with the 
electron tracking and detection efficiencies. 

Model dependence is estimated from the observed change in the moments when the 
B — > Dev and B — > D*ev decay shape parameters are varied according to their uncertainties. 
The uncertainty in secondary (B — > D — ► e) decays is derived from switching the branch- 
ing fraction corrections on and off. Similarly, the correction to the prompt B — > D^eu 
branching fraction in the Monte Carlo, incurs such a systematic uncertainty. 

The uncertainty due to mis-tagging in the B° and B + samples is derived from the mag- 
nitude of the combinatorial background. Uncertainty in the measured luminosity adds to 
the uncertainty on continuum electron yield. Additionally we account for the uncertainty 
associated to magnitude of the hadron fake contribution. The systematic due to the overall 
fit for the secondaries to the data is estimated by varying the lower p* e B bound of the fit. 

To estimate unfolding uncertainty we vary the effective rank parameter in the SVD al- 
gorithm. The uncertainty due to the b — > u subtraction, which occurs after unfolding, is 
evaluated by varying the normalisation of the De Fazio— Neubert inclusive spectrum. 

The total systematic error is obtained by adding each contribution in quadrature. It 
is important to note that the systematic errors are limited by Monte Carlo statistics. The 
contributions to the systematic error are summarised in Tables ITTT1 II VI IVl for the first, second 
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and third moments respectively. 

TABLE II: Measured moments, Mi, M2, Ms and the branching fraction for B — ► X c ef in the S 
meson rest frame for six cutoff energies E cut . The first error is the statistical, and the second error 
is the systematic. The moments are corrected for QED radiative effects. 



£ cu t[GeV] 


Mi[MeV] 


M 2 [10- 3 GeV 2 ] 


M 3 [10- 6 GeV 3 ] 


0.4 


1397.7 ± 5.1 ± 5.4 


172.8 ± 2.4 ± 2.2 


-22.4 ± 2.3 ± 0.7 


0.6 


1431.8 ± 4.8 ± 4.3 


148.2 ± 1.9 ± 1.2 


-11.6 ± 1.7 ± 0.6 


0.8 


1481.0 ± 4.4 ± 3.4 


119.9 ± 1.6 ± 0.9 


-3.9 ± 1.2 ± 0.6 


1.0 


1550.8 ± 4.0 ± 2.9 


89.0 ± 1.2 ± 0.4 


0.5 ± 0.8 ± 0.3 


1.2 


1631.6 ± 3.6 ± 2.2 


61.9 ± 0.9 ± 0.6 


1.9 ± 0.5 ± 0.2 


1.5 


1775.8 ± 3.0 ± 2.3 


29.4 ± 0.6 ± 0.3 


1.6 ± 0.2 ± 0.1 



TABLE III: Breakdown of the systematic errors for the first moment, Mi, for B — ► X c eu in the B 
meson rest frame for all 6 values of E cnt 





Mi 


Mi 


Mi 


Mi 


Mi 


Mi 




[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


E cut [GeV] 


0.4 


0.6 


0.8 


1.0 


1.2 


1.5 


electron identification 


0.69 


0.59 


0.57 


0.50 


0.44 


0.27 


detection efficiency 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


B — > D^eu form factors 


2.08 


1.53 


1.06 


0.89 


0.84 


1.07 


B -» D^eu Br 


4.06 


3.69 


3.04 


2.61 


2.00 


1.50 


B -> D 1 :*} -> e Br 

(s) 


0.21 


0.18 


0.13 


0.08 


0.05 


0.01 


continuum background 


0.35 


0.35 


0.30 


0.19 


0.08 


0.04 


combinatorial background 


0.01 


0.06 


0.06 


0.04 


0.03 


0.06 


hadron fakes 


0.70 


0.60 


0.45 


0.25 


0.08 


0.02 


secondaries 


1.99 


0.79 


0.67 


0.49 


0.01 


1.30 


unfolding 


1.63 


0.67 


0.32 


0.44 


0.31 


0.35 


b — > u subtraction 


0.17 


0.19 


0.21 


0.21 


0.21 


0.21 


total systematics 


5.38 


4.27 


3.37 


2.94 


2.20 


2.30 



A set of statistical correlation coefficients for all measured moments with E cnt ranging 
from 0.4 GeV to 1.5 GeV are given in Table IVT1 

SUMMARY 

We report a measurement of the electron energy spectrum of the inclusive decay B — > 
X c eu and its first, second and third moments for threshold energies from 0.4 GeV to 1.5 GeV. 
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TABLE IV: Breakdown of the systematic errors for the second moment, M 2 , for B — > X c eu in the 
B meson rest frame for all 6 values of E cvA 





M 2 


M 2 


M 2 


M 2 


M 2 


M 2 




[iu icrev j [iu 


(jrC\ J [1U 


(jCX J [1U 


(jrC\ J [1U 


-3r<a\r2~\ h n- 
ucv j [iu 


Lie v j 




n a 

U.4 


n p. 

U.D 


U.o 


1 n 

l.U 


1 9 


1.0 


electron identification 


0.31 


0.26 


0.2 


0.15 


0.11 


0.05 


detection efficiency 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


B — > D^eu form factors 


0.81 


0.43 


0.32 


0.29 


0.27 


0.16 


B -> D^eu Br 


0.18 


0.07 


0.22 


0.22 


0.22 


0.10 


B — > — > e Br 


0.05 


0.03 


0.02 


0.01 


0.00 


0.00 


pfrritmn 1 1 tti nn PKcrym inn 

VjUIlUlllU.U.1.11 UCH-'lVg ) l <J 1 11 ill 


0.12 


10 


06 


03 


02 


00 

\J . \J\J 


combinatorial background 


0.08 


0.06 


0.04 


0.03 


0.02 


0.01 


hadron fakes 


0.17 


0.12 


0.07 


0.03 


0.01 


0.00 


secondaries 


1.36 


0.78 


0.72 


0.13 


0.40 


0.09 


unfolding 


1.43 


0.78 


0.25 


0.16 


0.20 


0.21 


b — > u subtraction 


0.13 


0.11 


0.09 


0.08 


0.07 


0.06 


total systematics 


2.18 


1.19 


0.88 


0.42 


0.57 


0.31 



This set of moments, combined with the measurements of the semileptonic branching fraction 
and the moments of the hadronic mass distribution, will be used for the determination of 
the HQE parameters and of \V c b\. 
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FIG. 3: First, second and third electron energy moments, Mi, M2 and M3, as a function of cutoff 
energy E cnt . The errors shown are statistical and systematic. 
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